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Abstract 

The understanding of the spreading of hquids on sohd surfaces is 
an important chaUenge for contemporary physics. Today, the motion 
of the contact hne formed at the intersection of two immiscible fluids 
and a sohd is stiU subject to dispute. 

In this paper, a new picture of the dynamics of wetting is offered 
through an example of non-Newtonian slow liquid movements. The 
kinematics of liquids at the contact line and equations of motion are 
revisited. Adherence conditions are required except at the contact 
line. Consequently, for each fluid, the velocity field is multivalued at 
the contact line and generates an equivalent concept of line friction 
but stresses and viscous dissipation remain bounded. A Young-Dupre 
equation for the apparent dynamic contact angle between the interface 
and solid surface depending on the movements of the fluid near the 
contact line is proposed. 



Key words: contact angle, contact line, dynamic Young-Dupre equation, 
wetting 

1 Introduction 

The spreading of fluids on solid surfaces constitutes a significant field of 
research into the processes met in nature, biology and modern industry. In- 
terfacial phenomena relating to gas-liquid-solid systems take into account 
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contact angles and contact lines which are formed at the intersection of two 
immiscible fluids and a solid. The interaction between the three materi- 
als in the immediate vicinity of the contact line has great effects on the 
statics and dynamics of flows (Dussan, [6]). Many observations associated 
with the motion of two fluids in contact with a solid wall were performed 
(Bataille, [1]; Dussan and Davis, [5]; Dussan, Rame and Garoff, [7]; Pomeau, 
[24)). According to the advance or the recede of a fluid on a wall, we ob- 
serve the existence of an apparent dynamic contact angle when a contact 
line is moving. This angle, named after Young, depends on the celerity of 
the contact line, and the motion in the vicinity of the contact line does 
not seem to be influenced by the behaviour of the total flow (Bazhelakov 
and Chesters, [2]; Blake, Bracke and Shikhmurzaev, [3]). It is noteworthy 
that since Young's article on capillarity, [37], the understanding of these 
phenomena has remained incomplete. For example, it is well known that 
for Newtonian fluids the total dissipation and the interface curvature at 
the contact line are infinite (Huh and Scriven, [18]; Dussan and Davis, [5j; 
Pukhnachev and Solonnikov, [25]). In fact, fundamental questions remain 
unanswered. Among these are the following: 

What is the kinematics of the contact line? Can the fluid velocity fields be 
multivalued on this line? What is the work of the dissipative forces in its 
vicinity? Is there slip of the contact line on the solid wall? What is the 
connection between apparent and intrinsic contact angles? 
There are various ways to overcome these difficulties: to consider the slip 
length on the solid wall (Hocking, [15]; Shikhmurzaev, [32j). to consider one 
phase as a perfect fluid, the possibility of a thin film as a precursor film on 
a wall (de Gennes, [9]), the assumption of dynamic surface tension different 
from the static counterpart (Shikhmurzaev, [32]), the use of non-linear cap- 
illary theories such as Cahn and Hilliard's theory of capillarity (Seppecher, 
|29j ) , or the direct computation of flows by means of molecular models (Ko- 
plick, Banavar and Willemsen, [20]). All these attempts are not able to 
produce a complete satisfactory answer to the previous questions. 
It was noticed by using molecular methods that large amplitude shearing 
rates reveal a tendency to reorganize the liquid, to facilitate the flow and to 
reduce the viscosity. This suggests that in reality there may be rheological 
anomalies around the contact line (Heyes et al, [H]; Holian and Evans, }16j ; 
Ryckaert et al, [28]). 

For condensed matter and far from critical conditions, interfaces which are 
transition layers of the size of a few Angstroms between fluids or between 
a fluid and a solid can be modelled by surfaces endowed with a capillary 
energy (Rowlinson and Widom, [27J). Solid walls are rough on a molecular 
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or even microscopic scale. Moreover, the chemical inhomogeneity due to the 
nature of the solid or the presence of surfactants changes the surface tension 
in a drastic way. Nevertheless, roughness for example is taken into account 
by corrections of the measurement on a mean geometric surface (Wenzel 
equation in Cox, [Ij or Wolansky and Marmur, [36j). 

The motion of liquids in contact with a solid wall will be considered in this 
paper within the framework of continuum mechanics. Knowledge of the 
equations and boundary conditions which govern the movements of liquids 
in contact with solid walls and control of the contact line motion are the 
aim of our study. 

We propose a model of the dynamics of wetting for slow movements. To 
prove its accuracy, we are only considering partial wetting when the bal- 
ance contact angle is theoretically defined without ambiguity (de Gennes, 
Brochard-Wyart and Quere, |10j). The liquids are non-Newtonian; so the 
viscous stress tensor deviates from the Navier-Stokes model for large values 
of the strain rate tensor. For two-dimensional flows, and in the lubrica- 
tion approximation, the streamlines have an analytic representation and it 
is possible to obtain the flows near the contact line. Equations of motion, 
boundary conditions and some consequences on the contact angle behaviour 
are deduced. 

The notation is that of ordinary Cartesian tensor analysis (Serrin, [30]). 
In a fixed coordinate system, the components of a vector (covector) a are 
denoted by a*, (oj), where i = 1,2,3. In order to describe the fluid motion 
analytically, we refer to the coordinates x = {x^,x'^,x^) as the particle's po- 
sition (Eulerian variables). The corresponding reference position is denoted 
by X = (X^, X^, X^) (Lagrangian variables). The motion of a fluid is clas- 
sically represented by the transformation x = (f{t, X) or x* = 93* (t, X) . It is 
assumed that (f possesses an inverse X = x) and continuous derivatives 
up to the second order except at certain surfaces and curves. The vector 
V denotes the fluid velocity. The whole domain occupied by the fluid in 
Lagrangian variables is Dq and its boundary is the surface Sq. In Eulerian 
variables, the fluid occupies the volume Dt with boundary corresponding 
to the fixed regions Dq,T,q in the reference configuration. A moving curve 
Tf on Tit in the present configuration corresponds to the moving curve Toj on 
So in the reference configuration. The domains Dq, Dt,T,Q,Tit,TQt,^t must 
obviously be oriented differentiable manifolds. 

To each point of a unit normal vector n (n*), external to Dt, and a mean 
radius of curvature, Rm can be assigned. Furthermore Id is the identity 
tensor with components . Then Id — n (8) n (components 5^ —n^rii ) is the 
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projection operator onto the tangent plane of the surface S^; let t denote 
the unit tangent vector of Tt oriented; n' = n x t is the binormal vector to 
with respect to S^; it is a vector lying in the surface S^. 

2 General kinematics of a liquid at a contact line 

Following Dussan and Davis' experiments for contact line movements, [5], 
the usual stick-adhesive point of view of fluid adherence at a solid wall is 
disqualified in continuum mechanics. A liquid which does not slip on a 
solid surface does not preclude the possibility that at some instant a liquid 
material point may leave the surface. The no-slip condition is expressed as 
follows: 

The velocity of the liquid must equal the solid velocity at the surface. 




Figure 1: A liquid (in drop form) lies on a solid surface dS. The liquid 
La is bordered with a fluid Lb and a solid S; T,ti is the boundary between 
liquid La and solid S; is the boundary between fluid Lb and solid S, 
and consequently, S^i U SJ^ = dS; Et2 is the interface between liquid La 
and fluid Lb', ni and n2 are the unit normal vectors to Sfi and Sf2, exterior 
to the domain of liquid La and the domain of fluid Lb, respectively; the 
edge Tt (or contact line) is common to T,ti and 5]t2 and t is the unit tangent 
vector to F^ relative to ni; n'l and n'2 are the binormals to F^ relative to 
Eti and Sj2, respectively. 
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Let a liquid be in contact with (i) a solid body S on an imprint Sti of 
the boundary dS of S and (ii) an incompressible fluid Lb along an interface 
^t2 (fig- E])- Let, moreover, the mobile surfaces 11^ (a = 1,2) be described 
by the Cartesian equations fa(t,x) =0 (a = 1,2) and let the equations of 
the common curve Tt = Sti Pi Sf2 be given by 

/l(t,x) = /2(t,x) =0. 

For a geometric point M of Y,ta with velocity we obtain the kinematic 
relation 

in which the usual convention that over a doubly repeated index summation 
from 1 to 3 is understood. With the notations of fig. [H if we observe that 
n-ai = ^a{dfa/dx^) where is a suitable scalar, the celerity of the surface 
T,ta has the value Cq = riaiW^ = —Xa{dfa/dt). This celerity depends only 
on the coordinates (t, x) of M. 

Let the velocity of a point of Tf be denoted by W. We shall denote the unit 
tangent vector to Tt relative to ni simply by t and consequently, n'j^ = ni xt, 
= n2 X t (^^6 fig- 1)- l^fi^ velocity of the common line is then expressible 

as 

u = (Id - t t)W. 

It is orthogonal to Tt and its expression depends only on the coordinates 
{t, x) of the point on Tt but it is not necessarily tangential to dS. Then, u = 
/3in[ + (320.2, where Pi and P2 two scalars. Thus, along Tt, riaiW = Ca 
and ci = niiW^ = nuu^ = /32(ni, n2, t), C2 = n2iW'^ = n2iU^ = /3i(n2, ni,t). 
Consequently, 

C2 n'l - ci n'n 

u = =■ 

(ni,t,n2) 

in which (a, b, c) is the triple product of the three vectors a, b, c. Due to 
the definitions of ni and n2, we remark that (ni,t, n2) > 0. 
The kinematics of fluids in the vicinity of the contact line will be axiomatized 
as follows: 

T,ti is a part of the surface of the solid S. Liquid La adheres to dS in the 
sense of the no-slip condition previously proposed. St2 is a material surface 
of liquid La- 

At the contact line, the velocity u may have any direction between ni and 
n2, depending upon, how the contact line is approached within La- On the 
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solid surface Sfi of S it is, however, tangential to T,ti, so that 

U = — - = — i- = -U 111 

(ni,t,n2) (ni,t,n2) 

where u denotes the value of the contact line celerity in the direction liquid 
La to fluid Lb and V2 is the common velocity of the fluids on S(2- The 
contact line Tt is not a material line of La] its velocity is different from the 
velocities of the liquid on S and on S(2- 




Figure 2: Typical two-dimensional motion of fluids in contact on a solid 
surface with a stationary contact line. The wedges formed by Sji, 5]j2 and 
T,t2 , '^ti bound the fluids La and L b ■ The auxiliary surface S^s separates L b 
into two domains. A control surface together with S^i and S^^^ constitute 
the boundary of a compact domain Dt of the two fluids. For explanations, 
see main text. 

The motion of the particles of La on dS and Sj2 is comparable with that of 
an adhesive tape stuck on a wall, the other edge of the adhesive tape being 
mobile (fig. [2]): for u > (or C2 < 0) the particles of La belonging to T,t2 
are driven towards Tt and necessarily adhere to S along Tm- For u < (or 
C2 > 0) the result is reversed: the particles of La belonging to S^i reach Fj 
and are driven towards T,t2- 

In fig. [3] the motion of the fluids is sketched. The two manifolds S^i and St2 
constitute two sheets of the same material surface. The motion of the liquid 
La is represented by using a continuous mapping (p from a half reference 
space Dq{La) bounded by dSo, see fig. 3, onto the actual domain Dt{LA) 
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occupied by L^. The domain Dt{LA_) is included in the dihedral angle 
formed by S^i and Sj2- The contact line Tt is the image of the mobile curve 
Tot on dSo- Outside Tot, the mapping ip is C^-differentiable. 



Figure 3: In the reference configuration, the two sheets of the same material 
surface of fluid La are represented by a manifold OSq differentiable along 
Tot- Its image in present configuration, Dt is divided into two differentiable 
manifolds Sti,Sj2 forming a dihedral angle. The common edge Fj is the 
image of a moving curve Fot in dSo. The triad no,to,nQ, in the reference 
configuration is transformed by the mapping (f into the triad ni,t,n'^ or 
the triad n2,t,n2 depending on whether its image is associated with the 
manifold S^i or the manifold T,t2- For other notations see main text. 

A second fluid Lb occupies the supplemental dihedral angle (Sji, ^t2)- The 
conditions of motion are the opposite to those of L^. The material surface 
T,t2 is the common interface between La and Lb- Provided the fluids are 
not inviscid, the velocities of the fluids La and Lb are equal along Y^t2- 
Moreover, for liquid La, if it > 0, the particles of La are driven towards Fj 
and adhere to S along Y^ti- The particles of Lb are also driven towards Fj; 
if they adhere to S along S^^, the contact line goes through the two fluids 
La and Lb along the solid wall dS, (if n < 0, a change in the time direction 
along the trajectories leads to analogous consequences). This is in direct 
conflict with the fact that Fj belongs to the interface T,t2 separating La and 
Lb- To remove this contradiction, it is possible to separate the fluid Lb iti 
two parts with a material surface T,t3 (see figs. [2]and[3|). S^]^ and T,t3 are 
the two sheets of the same material surface for a domain D[[Lb) of the fluid 
Lb within the wedge formed by the dihedral angle (J^'n, ^^a)- The sheets 




n 



reference configuration 



present configuration 
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T,t2 and T,t3 constitute two parts of the same material surface for a domain 
D'1{Lb) of the fluid Lb within the wedge formed by the dihedral angle 
(St2,St3). The two domains D[{Lb) and D'/{Lb) with common material 
surface Sj3 - across which velocity is continuous - constitute two independent 
fluid domains which do not mix. For the domain D^{Lb), the conditions in 
the vicinity of the contact system are similar to those of liquid L^. Velocities 
which are discontinuous and multi- valued on Tt, are compatible with the 
movements of fluids and Lb within the domains Dt{Lji), D[{Lb) and 
D'I{Lb). 

3 Equations of motion and boundary conditions 
revisited 

The fundamental law of dynamics is expressed in the form of the Lagrange- 
d'Alembert principle of virtual work applied to any compact domain of the 
two fluids La and Lb- 

Any compact domain Dt made up of the two fluids La and is bounded by 
Tjt- The boundary Et is constituted of S^i, S^j^ and a complementary surface 
Sj which is not in contact with the solid surface dS: T^t = '^ti U U Sj 
(fig. 2). In a Galilean frame, the virtual work due to the forces applied to 
La and Lb (including inertial forces but without forces due to capillarity) 
is in the general form 

/ [{cl,,-pai)C + {p5{-Qi)C,j]dv+ f nCda. (1) 

JDt Jj^t 

Here, dv, da (and later dl) are the volume, area (and later line) increments, 
C denotes any virtual displacement field, the volumetric forces, p the 
density, a the acceleration vector, Q the viscous stress tensor and p the 
pressure. Moreover, the stress vector P describes the action of the external 
media on Sj. A contribution along Tt is not accounted for. 

In continuum mechanics, fiuid-fiuid and fluid-solid interfaces are differen- 
tiable manifolds endowed with surface energie^. We denote aAB , (^AS and 
(Tb5, the surface energies of interfaces liquid L^-fiuid Lb, liquid L^-solid 

^In statistical physics, fluid interfaces are transition layers of molecular size. They are 
modelled in continuum mechanics with regular surfaces (Rowlinson and Widom, [27] )■ On 
a molecular scale, a solid wall is rough; but in continuum mechanics, when the scale of 
the roughness is vanishingly small relative to the size of the solid wall, the solid wall and 
the fluid-solid surface energy are modelled with a differentiable average surface, fiat on a 
microscopic scale and a corrected surface energy (Wolansky and Marmur, [36)). 
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S and fluid L^-solid S, respectively. It is usual to define a measure of en- 
ergy on interfaces denoted by a da, where a stands for ctaBjO'as or ass 
following the interfaces between the fluids and the solid. The total energy 
of capillarity of the interfaces Sji , S^^^ and T,t2 is 

E = I a da. 

For any virtual displacement field, the variation of E is (Gouin and Kosihski, 

= j '^^ - + {^l - "'■''^*) ^j) C' da 

+ j {{'yAS-crBs)nu + (TABn2?j C dl, 

where the scalar 6a is the variation of the surface energy a associated with 
the displacement vector n (n*) and scalar Rm stand for the unit normal 
vector and the mean radius of curvature to Eji, or St2) respectively. 
The surface energy uab between the liquid La and the fluid Lb is pos- 
itive and constant (Rowlinson and Widom, [27]). Generally, a fluid-solid 
surface energy depends on the fluid which is in contact with the solid, the 
geometrical and physico-chemical properties of the solid, the microscopic 
asperities or the presence of a surfactant. The simplest case occurs when 
the surface energy is defined as a function of the position on the surface 

(xG St ^CT(x,t)). 

Hereafter considering such a case, the virtual work due to the forces of 
capillarity applied to La and Lb is simply 

/ 4~ "-i C da - I ({aAS - (^Bs) n'li + gab n'2i) C dl, (2) 

and relations (1), (2) lead, after execution of the variations and performing 
integration by parts in several volume and surface terms, to the expression, 
denoted by bT , of the virtual work by forces applied to the domain Dt 

ST= [ [i^i-pai-p,^ + Ql) C]dv 

J Dt 

+ f (^Pi+pui- Q{ nj) C da + I (Pi + + p) nu - Q{ J C da 



+ {{ +PA- Pb) n2i - (Qai - Qbi) ^2j^ C da 
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-j {{aAS - (^Bs)n'ii+ (TAB n'2i^ C dl. (3) 

Unit normal vectors rii and n2 are exterior to the domain of liquid La and 
fluid Lbi respectively; as is called a as or ubs depending upon which fluid 
is in contact with 95, and Ql is called or Q'^j, respectively depending 
upon which fluid is in contact with Sj and dS. We emphasize that it is not 
necessary for Z^t, and Fj to be material. Finally, virtual displacements 
are tangential to the solid surface dS (on dS, nuC = 0). 

The expression of the Lagrange-d'Alembert principle is (Germain, ^llj): 

For any C, such that on dS, nijC* = 0, then 6T = 0. 

We emphasize that this principle is not associated with a variational ap- 
proach and there is no variational principle in it: 6T is not the Frechet 
derivative of a functional (Gurtin, [13j). Only for equilibrium, and due to 
the fact that the viscous stress tensor is null, the minimization of energy 
(this is a variational principle) coincides with this approach. Such a method 
is relevant to the theory of distributions where ^ are vector fields of class 
C°° with compact support (Schwartz, [31j). 

The equations of motion and natural boundary conditions that emerge from 
it are as follows: 

Equations of motion 

pa^+p,i = (p. + Qlj. (4) 
Conditions on the liquid La - fluid Lb interface St2 

^^^n2j = {Qai - QsM^j + (PB - PA)n2i, (5) 
which is the dynamic form of the Laplace equation. 

Conditions on the boundary 

Pi = -pni + Qluj, (6) 
which is the classical expression of the balance of stresses for viscous fluids. 

Conditions on the surface dS = Y^ti U 
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Expression (3) of the virtual work and the Lagrange-d'Alembert principle 
imply: For any C such that on dS, nuC = 0, 

/ + i^+P) - Qi ^ij) C da = 0. 
Jas \ J 

Consequently, there exists a scalar field x of Lagrange multipliers defined 
on dS such that (Kolmogorov and Fomin, [19j) 

Pi = Qi nij - {'^^ + p) nii + X nii. (7) 

Generally P is not coUinear to ni and x is ^in additional unknown scalar. 

Conditions on the contact line Fj 

Due to the condition on dS, riuC = 0, a virtual displacement is expressed 
at any point of the contact line F^ in the form 

C = Kt + vn'i, (8) 

where the two scalar fields k and v are defined on F^. For any field ^ in 
the form (8), the contribution of /p^ (^{o'as — ^^Bs) ^'u + '^AB C dl in 
relation (3) yields 

j (^{(^AS - (^Bs) + GAB {Kf + v n'l ) dl = (9) 

In the general case, since n'l = ni x t and n'2 = n2 x t, expression (9) 
implies 

(^AS-cTBS + a AB cos 9i =0, (10) 

where 9i is the angle between n'l and n'2. This angle named intrinsic contact 
angle in the literature (Wolansky and Marmur, |36j). is the angle in a plane 
that is normal to dS and F^ between tangents in O parallel to and dS 
(see fig. m). 

At point O of the contact line Fj, let us consider the section of the La-Lb 
interface in the plane erected by n'l and n'2. Let the curvature of the planar 
section Ct2 of St2 be . For a two-dimensional flow, the mean curvature 
of the surface St2 is 2i?~^ = R~^. At a generic point of Ct2-, the angle 
between n'l and the tangent to Ct2 is denoted by 6. This angle depends on 
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as 



o 




Figure 4: Cross section Ct2 of the fluid-fluid interface and the soHd wall in 
the plane (O, n'l, n'2). The contact line is reduced in this figure to the point 
O. The intrinsic angle 6i is the angle between n'l and n'2. At any generic 
point of Ct2, the angle between n'l and the tangent to Ct2 is denoted by 
9; thus, at the points M, P and A, the angle 6 is denoted by 9p and 
9a, (see main text in section 3). The intersections of the tangent lines to 
Ct2 in M and P with the axis On'i are denoted by and Ip. The arcs 
of circles C{Im) and C{Ip) of centers Im and Ip and radius ImM and IpP 
intersect the axis On'i at points Jm and Jp, respectively. A point Q of 
C{Ip) is represented by (^,0;) in the polar coordinate system of pole Ip and 
polar axis On'i (see main text in subsection 5.2). 

the choice of the point and on the fluid flow. Since the surface energy aAB 
between two fluids is constant and R = dl/d9, we obtain, 




O'AB 

R 



sin 9 dl = a AB {cos 9i — cos9p), 



(11) 



where 9p is the value of 9 at the point P. Relation (5) yields 



O'AB 

R 



{Qm - <3si)"'2"2j +PB- PA, 
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and, consequently, relations (10), (11) yield d 

cfAS - CFBS + CFAB COS + J ((^ii " Qsi) "'2"-2i + Pb - PA^ sin e dl = 0. 

(12) 

4 A creeping flow example of non- Newtonian fluids 
4.1 The Huh and Scriven model revisited 

To understand more precisely the behaviour of a liquid near a moving solid- 
liquid-fluid contact line, we reconsider the situation of two-dimensional flows 
proposed by Huh and Scriven, [18]. Let us recall the main results of their 
article (see fig. [2]): 

A flat solid surface in translation at a steady velocity U is inclined from a 
flat interface between a liquid La and a fluid Lb (here, the angle 9 of fig. 
4 is constant, independent of the generic point of Ct2)- The contact line 
velocity with respect to the solid is — U; thus in the notations of section 2, 
we obtain U = n In a two-dimensional situation of the plane (O, n'^^, rig), 
it is convenient to take the contact line intersection point O as the origin of a 
polar coordinate system (r, (p) and O n[ as the reference polar axis. The two 
bulks are incompressible Newtonian fiuids. In term of the stream function 
ip{r,ip) for two-dimensional steady flows, the velocity is 

, 1 dip dip 

V = VrGr + U<pe<p, With Vr = -"T^, = T^, 

r Oif or 

and (0,6^,6^) as the mobile polar frame. In the creeping flow approx- 
imation of a viscous fluid, Eq. (4) leads to the linearized Navier-Stokes 
equation and consequently to the biharmonic equation V^ip = , where 
and are respectively the Laplacian and bi-Laplacian operators (Moffat, 
|23j : Bataille, [1]). A solution is in the form: 

ip{r,ip) = rf{ip) 

which leads to the ordinary differential equation 

/(^)+2r(v.)+/(^^)((^) = o, 

^ Relation (12) can be proved directly by using the projection on n'l of the balance of 
forces applied to a liquid La - fluid Lb domain containing the fluid interface Et2. 
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and thus to the general solution 

ip{r, if) = r{a sin ip + h cos ip + c sin ip + d ip cos (p), (13)i 
which holds for either fluid. 

The boundary conditions at the solid wall and the liquid-fluid interface are: 
(i) a vanishing normal component of the velocity at the solid surface and 
interface, 

(m) continuity of the velocity at the interface, 
[in) continuity of the tangential stress at the interface, 
[iv) non tangential relative motion of the fluids at the solid surface except 
at the contact line. 

These eight linear conditions yield the values of coefficients o, 6, c, d for the 
two fluids La and Lb- If the dynamic viscosity coefficients are identical, the 
eight integration constants are: 



OA 


= -ueD{e)[ 


TT — 6 + sin 6 cos 9 ] , 


as 


= uD{e)[ (vr 


- 6l)(sin 61 cos 61 - 61) + vr 6* sin^ 


hA 


= 0, 




bB 


= uD{e)[ vr i 


^ sin 6 cos — TT sin^ ], 


< 

CA 


= uD{e) (vr- 


- 0) sin^ 0, 


CB 


= -ud[9) e 


sin^ 0, 


dA 


= uD{6)[ (vr 


- 0) sin 6* cos 6* + sin^ ], 


. d-B 


= uD{e)[ -6 


' sin cos + sin^ ] 



with D{0) = [ 61 (vr - 61) - (vr - 261) sin cos - sin^ 
No difficulty should arise, in principle, in the determination of ua, ■ ■ ■ ,dB 
for two fluids with distinct dynamic viscosity coefficients. The form of the 
streamlines as obtained by Huh and Scriven are sketched in fig. [2l and the 
motion of the contact line fits perfectly with the kinematics as outlined in 
section 2. Furthermore, if /i^ = fJ-B = ^ values of the dynamic contact 

angle 0, the viscous stress components are 

Trip = - -(/ + /") = - 2 -(c COS ip - d sin ip), Trr = T^^ = 

r r 
and the pressure field is given by 

P = Po + -{f' + f") = - 2 -(c sin +d cos p), 
r r 

where po is the hydrostatic pressure (in both formulae phase subscripts have 
been omitted). 
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As proved by Huh and Scriven, the dissipation in any domain Dt of the fluids 
containing the contact hne Tt, ^ (/ + f")'^ /r"^ dv^, and the total 

traction exerted on the soHd surface by the fluid interface are logarithmicaUy 
inflnite. Moreover the normal stress across the fluid interface varies as r~^; 
furthermore, the stress jump should be balanced by the Laplace interfacial 
tension uab and the curvature does increase indefinitely at the 
contact line. These are all non-integrable singularities. 

4.2 A model of non-Newtonian fluid near the contact line 

To avoid the previous paradox of an infinite dissipative function at the con- 
tact line, we consider non-Newtonian incompressible fluids with a convenient 
behaviour of the viscous stress tensor. It is experimentally known that the 
dynamic viscosity /i of polymeric liquids depends on the shear rate e. The 
behaviour prevailing in such situations is not well understood. In a wide va- 
riety of technological applications, liquids are subjected to large shear strain 
forces. A molecular dynamic investigation of liquids subjected to large shear 
strain rates has been performed by Heyes et al, [M] . The shearing action has 
been found to change the liquid structure and reveals a tendency to reduce 
the shear viscosity (Ryckaert et al, [28\). 

In the literature, some empirical formulas for the viscosity obtained by 
means of a weighted least-squared adjustment were proposed. For exam- 
ple, fj,{e) ~ ^(0) — c was suggested as possible form for the viscosity 
as a function of e (Heyes et al p^). Holian and Evans, [16], proposed a 
representation in the form 

/i(e) ~ fi{0) - c^e. (14) 

Data for the viscosity of an atomic fluid generated by nonequilibrium molec- 
ular dynamic were performed by Ryckaert et al, [28]. They indicated that 
for shear rates below lO^^s"^, /i(e) does not differ significantly from /i(0) 
but that these previous laws are not extendable when e tends to infinity. 
We propose a model where the viscous stress tensor Q is a function of the 
strain rate tensor A = (1/2) (VV -|- (VV)*). For moderate values of A the 
fiuid is Newtonian and the function is linear. The function deviates from this 
classical behaviour for large values of A. For an isotropic stress tensor of 
two-dimensional flow, the Rivlin-Ericksen representation theorem (Truesdell 
and Noll, [34]) leads to a viscous stress tensor in the form Q = A I -|- 2 fiA. 
but A and /i are non-constant functions of invariants of A, and /iA is a 
functional of A where tends to //g (a*o being constant) when A tends to 
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zero. 

We propose to use a convenient representation of the /x-behaviour in the 
form ^ 

/X = (1 - e-^) /xo with 7 = (15) 

where || A|| = is the norm of the strain rate tensor, tq is a characteris- 
tic time of the fluid and ^ is a small parameter (0 < ^ <C 1). Then, for very 
high shear rates /i behaves as a step function as expected in Ryckaert et al. 
To fit with relation (14) when e is close to lO^^s"^, we choose ^ = 0.2 
and To = 10~-'^^s but many other values can be considered and results of the 
literature are disparate. 

In the following, we take ^ = 0.1 and tq = lO^^^s ; then for 7 = 7q = 4.6, 
we obtain /x = 0.99 //q- We call this 70-value the cut-off coefficient. For 
7 = 7q, ;U ~ /Uq, and the fluid may be considered as Newtonian. 

For the Huh and Scriven model of two-dimensional incompressible flows, 
ll^ll = e = |/ + /"I /(2 r). When u = 1 mm. s~^, due to the 79— value, 
considerable variations of /it occur from the contact line to a distance of 20 
to 30 Angstroms. The same holds true on the solid wall when G [5 , 175 ] . 
Outside these distances from the contact line, the fluids can be considered 
as Newtonian. Then, tends to zero for very large values of the shear rate, 
and Q is a function of A which tends to infinity with A, but weaker than 
a linear function. The total stress tensor of a fluid is always of the form 
— p Id + 2 yuA where p = 11 — A (here 11 notes the hydrostatic pressure). 
For steady flows, the equation of motion is 

p grad (^ V^) + protV xV + grad p = n V^V -|- 2 A.grad n + pg (16) 

where g denotes the acceleration due to gravity. For a stream function 
'0(r, (f) = rf{(f), we obtain 

= ((/' + /-)e, + (/ + f)e^) , 

2 A.grad /x = ^ Mo ^ ((/' + f"')er - (/ + f")e<p) , 
p rot V X V = - ^ (/ + /") [fer + f'e^) . 

The inequalities < 7 e~'^ < 1 — e^^ , < ^ <C 1 and the fact that 

near the contact line 1 <^ --^^ — -—r yields the approximate form of Eq. 

rp if + f") 
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(16) ^ 

p grad (2^^^) + P = /^V^V + p g, 

which imphes 

/X rot(V2v) + grad n x V^V = (17) 

with 

^l rot(V2v) = - ^ (1 - e-^) (/ + 2f" + /(^^)) k 

and 

grad ^ X V V = e /xo [f + f + ^ ^ j„ j k 

where k denotes the normal vector to the plane of the flows. 
From ^ 7 <C 1 — e~'^, we deduce again the Huh and Scriven ap- 
proximation for the stream function in the form tp{r,ip) = rf{ip), with 
/((^) + 2/"(v9) + /(^^)((^) = 0. 

On the solid wall, adherence conditions are required. This assumption is 
in agreement with molecular dynamics of fluid flows at solid surfaces: The 
non-slip boundary condition appears to be a natural property of a dense liquid 
interacting with a solid wall with molecular structure and long range force in- 
teractions (Koplick, Banavar and Willemsen, t20j). For the non-Newtonian 
model and for the creeping flow approximation, the general solution (13)i 
holds true for either fluid. Furthermore, the boundary conditions at the solid 
wall and at the liquid-fluid interface are the condition (i) — (iv) of subsection 
4.1. Consequently, for the non-Newtonian model and for the creeping flow 
approximation the trajectories and the velocities near the contact line are 
identical to those of fluids with constant viscosity /Uq in the Huh and Scriven 
model. 



The dissipation H in the domain Vi = {r £ [0,1], ip G [0,6*], z G [0, L]} 
(where z denotes the contact line coordinate) is 

H = / tr(/.A2)d^ = L^ f f (/ + f"? dr d^. 

Jvi ^ Jo Jo r 



From the inequality 0<1 — e '^<7, we deduce 

^ ^ Jo Jo tI ^ 2i~€ tI Jo 
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which proves that the dissipation is finite at the contact Une, since / given 
by (13)i, (13)2 is bounded. 

Other calculations yield a bounded total force exerted on the solid surface by 
the fluid-fluid interface near the contact line, but another problem associated 
with the fluid-fluid interface curvature still remains unresolved. 



5 Study of a curved interface in the vicinity of the 
contact hne 

For a Newtonian fluid in two dimensional flows, the fluid interface curvature 
should increase rapidly as the contact line is approached. This result is 
in direct conflict with the hypothesis of section 4 that the fluid interface 
is perfectly flat. Indeed, Huh and Scriven, [H], pointed out when water 
at moderate dynamic contact angle wets a surface at 6 mm.min~^ , the local 
radius of curvature would have to be about 10^ time greater than the distance 
to the contact line and the curvature would be imperceptible by optical means. 
Nevertheless in such a case, the intrinsic angle 6i at the contact line may 
strongly deviate from the angle 6p, (see fig. 4), which is observed at a point 
P near, but not at, the contact line. 

Let us consider results presented in subsection 4.2; the equation of motion 
(4), boundary condition (5), and calculations of subsection 4.2 yield the 
pressure field values for the fluids La and Lb, 

PB-PA = ^ iif'B + Ib) - if A + /a)) ^ — 2vrsin 9 D{e). (18) 



r 



As done by Huh and Scriven, we notice that 

(^AB 



and consequently, 



PB-PA ^ 
r fi u 



2TTsm9D{e). (19) 

R GAB 

For partial wetting, it is easy to compute the value of 2tt sm.9 D{9) numer- 
ically; we obtain 

-<9<— 4<2TTsm9 D(e) <15. 

6 6 

When the capillary number Cq = /ig \u\/aAB is sufficiently small, (in ex- 
periments, Ca is often smaller than 10~^), we deduce \r/R\ ^ 1; it is all 
the more true, for a non-Newtonian fluid given by the representation (15), 
where ^ tends to zero at the contact line. 
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5.1 Two-dimensional steady flows near the contact line 

The conditions and the notations are given in subsection 4.2, but the fluid- 
fluid interface is curved. The cross section of the fluid-fluid interface and 

the solid wall is presented on fig. 4. We consider the domain occupied by 
the two fluids in the immediate vicinity of the contact line and we assume 
that, along Ct2 

T T do 

lim — = Hm — - = 0. (20) 

r^O R r^Q dl 

Thus, in the immediate vicinity of the contact line, we obtain \r/R\ <C 1. 
On fig. 4, at the generic point P of Ct2, the intersection of the tangent 
line with the axis On'^ is denoted by Ip. To each point Ip corresponds, 
in the fluid domains, the arc of a circle, denoted by C{Ip), with center 
Ip and radius IpP. To a point Q of C{Ip) corresponds the polar coordi- 
nates, {g,Ld), associated with the pole Ip and the mobile frame (e^,etj). 
The polar coordinates of P are {g, Op) and will be denoted simply by [g, 9). 
Let us denote by y = gsinO, the distance from the point P to the solid 
wall. We deduce, dy = dg sin 9 + g cos 9 d9 and due to the diff'erential 
relation, dy = dlsm9, we obtain, dg/dl = 1 — g cotg6' d9/dl. When 
9 7^ or TT, limr^o g cotg 9 d9/dl = and consequently, when r — 0, 
dg ^ dl. Let M be another point of Ct2', we denote Jp and Jm the intersec- 
tions of C{Ip) and C{Im) with the axis On'^; when M ^ P, ||JmJp|| ~ dg, 
||lmlp|| ~ g d9/sm9 and relation (20) implies, limr^o(||ImIp||/|| JmJpll) = 0. 
Thus, in the immediate vicinity of the point O, HI^IpH <C ||JmJp|| and the 
two arcs of circles C{Im) and C{Ip) are distinct. In the following, we prove 
that, near the point O, a point Q of the fluid domains is represented by the 
orthogonal coordinate system {g,io). The equation of the curve Ct2 can be 
written in the form 9 = '&{g). 

As in section 4 for two-dimensional steady flow, the stream function ^[g,uj) 
verifies V^^* = 0. We look for a stream function in the form 

^{g,Lo) = gh{u,9), 

where 9 = '&{g), and such that the partial derivatives of h with respect to u 
and 9 are bounded. But, 

dOQ = dOIp + dIpQ = ( n[ + eg] dg + g dio e^,. 

\ sm 9 dg J 

Since lim^^o g d9/dg = and sin 9 ^ 0, we obtain 

dOQ = dg {eg + Oi{g)) + g duj e^^, with lim Oi(^) = 0, 
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and near the point O 

dOQ = dg eg + g dio e^^. 

Moreover, 

, ^ dh de, dh 

grad* = (/. + ,--) e, + -e.. 

and 

dh 

grad * = (/i + 02(£>)) + — e^j, with hm 02(£>) = 0. 
The velocity is 

dh 

and near the contact line 



V = ■Ug e^ + t^a; ©0,+ 03(^1) with, Vq = — — , Vu, = h and, lini 03(^1) = 0, 



In the following, On{g) with n ^ N, denotes a smooth function of the order 
of g such that lim^^o On{g) = 0. Similarly, 

V^* = div grad = - h+— —+- 7^3+- ^ = - h+- 05 ^ ^ 

g ou dg g ouj^ g Q Q c''^ Q 

^/r,2T^ ( 1 /, S^^N 1 ^dh d^h ^dO 1 , , 
grad(V2*) = (-^ + ^) + - + + ^ 

and finally, 



1 92^, 1 ,dh d^h ,de 1 ,^2/1 d^h, 1 , , 
g^ ^ dbJ^' g^ ^89 duj^dO' dg ^ ^duj^ doj"^^ g^ 

1 d'^h d^h, 1 , , 

= ?^^ + 'a:;^ + a7) + ?"«^^)' 

The principal part of the Laurent expansion in g of V^^* leads to the partial 
derivative equation 

d^h d^h ^ 

and the general solution for the principal part of ^ is 

'^{g,u}) = g {a sin lo + b cos lo + c lo sin to + d lo cos lo), (21) 
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where a, b, c, d are functions of 0. When = fj,^, the boundary conditions 
at the sohd wah and the conditions {i)-{iv) presented in section 4 yield the 
values of coefficients a,b,c,d for the two fluids and Lb- These values 
are given by relations (13)2, but here, 9 is not constant. The viscous stress 
components are 

Tguj = [h + —t) = - 2 - (c cos w - d sm uj), Tgg = t^jlu = 0. 

Q OLO'^ Q 

As in section 4, the pressure field is given by 

p = Po — 2 — (c sin CO + d cos w). 


In partial wetting, the stream function (21) together with relations (13)2 
show that the partial derivatives of h with respect to uj and 9 are bounded 
along Ct2- Along Ct2, Pa — Pb = (^ab/R and consequently, 

^ = J^2nsm9D{9). (22) 
dl a AB 

Furthermore, ||A|| = {1/(2^))} | h + d'^h/ doj'^\^=e- Thus, 

I T/ 1 1^1 I I 

II A|| = — \cA cos 9 — dA sin 9\ = — \cb cos 9 — ds sin 9\ = — sin^ 9 D{9). 
Q Q Q 

Let us note that, in the frame (O, n'^^, ni), relations x = q cos 9, y = q s'm9, 
dy = dl sm9, 9i = Arccos{((7B5 — (yAs)/(^AB} together with Eq. (22) and 
relation (15), allow us to obtain the parametric representation of Ct2 near 
the solid wall. 

Eq. (22) allows us to verify that limj._^o rd9/dl = and thus, the choice 
of the stream function ^ in the form (21) together with relations (13)2 is 
justified in the vicinity of the point O. Let us note that when it > (resp. 
ti < 0), is an increasing (resp. decreasing) function of the distance ?/ of a 
point of Ct2 to the solid wall. We do notice also that, whereas the curvature 
of Ct2 tends to infinity when the point O is approached, the stream function 
^ has the same form as the stream function proposed in subsection 4.2 
for a plane interface (a good example of such a curve is given, near x = 0, 

3 

hy y = |x|2 (Voinov, [55])). 

5.2 Apparent dynamic contact angle and line friction 

Let A be the point of Ct2 associated with the cut-off coefficient value 79, 
defined in subsection 4.2: the point A is at the border between the Newto- 
nian and the non-Newtonian domains of the fiuid fiows. We call apparent 
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dynamic contact angle 6a, the value of 6 associated with the point A (see 
fig. 4). Along the fluid-fluid interface, condition (Q'^j — Q^i) ^2"'2j = 0, 
together with relation (12) and pb — PA = (1/f?) fJ-u 27r sinO D{6), imply 

CTAS - (TBS + (TAB cos 9a + yu = 0, (23) 

with ^ 

u = I ^2tt sin^e D{9)dl. (24) 
Jo Q 

Relation (23) is a form of Young-Dupre dynamic relation for the apparent 
dynamic contact angle. We call z^, the line friction. It is easy to verify that 
the scalar is positive and of the same physical dimension as a dynamic 
viscosity. This result corresponds to the assumption in the article of Stokes 
et al, [33], in which they say that there is an additional viscous force on a 
moving contact line. Other expressions for the line friction have also been 
proposed (an attempt is done by a thermodynamic point of view in Fan, 
Gao and Huang, [8]). 

In the case of equilibrium, relation (23) yields the static Young-Dupre rela- 
tion (Levitch, ^) 

<TAS - (TBS + (TAB cos 6e = 0, 

in which 6e is the balance Young angle and 6^ = 9i = 6a. For any value of the 
contact line celerity, relations (23) and (24) yield, implicitly, the apparent 
dynamic contact angle da- With the formula (23), a simple explanation of 
a well-known experimental result (Dussan, [6]), may also be corroborated: 
with the advance of the contact line, u is positive and the apparent dynamic 
contact angle 6a is larger than the equilibrium angle 6e. This result is 
reversed when u is negative. 

5.3 Numerical investigations of the apparent dynamic con- 
tact angle and the line friction 

Hoffman, [17J, Legait and Sourieau, [21], Rame and Garoff, [26], and many 
other authors experimentally observe that, near the contact line, for slow 
motions, the apparent dynamic contact angle seems independent of the mi- 
croscopic distance to the solid surface. Let us verify numerically this obser- 
vation. 

Using the relations y = g sin 6 and dy = dl sin 9, relation (22) implies 

dy ^ (TAB 
y jj, u2tt sin 9 D{6) 
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Let us consider a point of Ct2 in the Newtonian domain of the fluid flows. 
Then, 

In(^) = r de, (25) 

Va Joa /^o 27r sm 6 D{0) 

where ya and y„ denote the distances of points A and N to the solid wall. 
In partial wetting, when 7r/6 < 9 < 5tt/6, then < 27r sin 6 D{9) < 15 and 
if -u > 0, 

O<0„-0, <15C, In(^). 

Va 

If we consider the case when y„ = 10^ ya, a crude approximation yields 
Sn — 9a < 140 Ca and thus 6n — 0a tends to zero with Ca- For example, 
when Ca = 10~^, we obtain 9n — Oa < 0.014 radian, (i.e. 0.8 degree), and 
the apparent dynamic contact angle seems independent of the distance of 
the point to the solid wall: in the lubrication approximation for two di- 
mensional flows, Eq. (23) expresses the behaviour of the apparent dynamic 
contact angle independently of any microscopic distance to the contact line. 
This result is in accordance with Seppecher's calculations [29] . 

Let us estimate an order of magnitude of the line friction. Along Ct2, Eq- 
(21) implies that, for each fluid, \ h + d'^h/duj'^\^^g = \u\2sm^e D{e). When 
^ = 0.1, a numerical computation yields 

-<^<— ^ 1.015 < /i + — K < 1.032 ki^. 

Taking into account the norm || A|| of the strain rate tensor along Ct2, rela- 
tion (15) yields a value of 7 such that 

^ 1.02 Vron/ ' 

and d'j/'j = ^ dg/g. Then, relation (24) allows us to obtain the value of the 
line friction 

u = — — / sm'^6' D(9) dj. 

Jo 7 

Since D{9) > 0, we obtain 



no 1 _ g-T 

sm^9rD(9r) / d-f, with 9r£[9i,9a 

Jo 7 
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in which Or is a convenient angle. Due to / (1/7)(1 — e^'^)d'y ~ 1, we 

Jo 



obtain 



^""^^ singer Di 



In partial wetting, a numerical computation implies 

^<^r<-^ ^ 0.68 < sin^^^ D{er) < 1.17, 
and consequently 

42 HQ<iy < 73 fi^. (26) 

Eq. (10) and Eq. (23) imply | cos 9a — cos 9i\ < i'\u\/aAB- Taking into 
account inequalities (26), we obtain | cos 9a — cos 9i\ < 73 Ca- In the partial 
wetting case, for 9 G (7r/6, 57r/6), then | sin ( {9a — ^i)/2) \ <73Ca and the 
derivative of sin^ 9 D{9) belongs to (—1, 1); we can deduce | sin^ 9^ D{0,) — 
siv? 9a D{9a)\ < 2 Arcsin(73 Ca). This crude approximation proves that for 
Ca sufficiently small, sin^ 6*^ D{9r) is close to sin^ 6'a D{9a) (for example, 
when Ca = 10"^, we obtain | sin^ 9r D{9r) - sin^ 9a D{9a)\ < 0.015). Re- 
placing D{9a) by its explicit expression sets the approximative relation of 
the line friction : 



e V e„ (vr - e„) - (tt - 2 9a) sin 9a cos 9a - sin2 ^ >' 



6 Concluding remarks 

In this paper, the dynamical problem of the contact of two non-Newtonian 
viscous fluids La and Lb with a solid was analyzed. Except at the contact 
line, these fluids were assumed to adhere to the solid and to each other. 
Then, the principle of virtual work allows us to obtain the governing equa- 
tions and boundary conditions. It was shown that the equations of motions 
and boundary conditions lead to streamlines near the contact line similar to 
those of a Newtonian fluid endowed with a dynamic viscosity which is, when 
the strain rate tensor tends to zero, the limit of the dynamic viscosity of the 
non-Newtonian fluid. The analyze of the stress tensor and the dissipative 
function near the contact line leads to several remarks and conclusions: 

For dissipative movements, a viscous stress tensor was added to the 
pressure term and in the expression of virtual work, the dissipative terms 
were distributed within the volume as Ql ■ and on the surfaces as —Qluaj- 
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In magnitude, the surface tension is comparable with the bulk stress despite 
the fact that the thickness of the intcrfacial layer, where the surface tension 
acts, is negligible compared to the characteristic length scale in the bulk: 
although the interfacial layer is very thin, the intermolecular forces which 
act on it and give rise to the surface tension are very strong so that the 
result is finite. Huge variations of the fluid velocity appear at the three- 
phase contact line. The only physical factor, which achieves to magnify its 
role, does not come from intermolecular forces but from the discontinuity of 
the velocity at the contact line and consequently from the viscosity along 
the contact line when the thickness of the contact line region tends to zero. 
In the immediate vicinity of the contact line, the viscous stresses yield a 
friction force which acts on the fluid-fluid interface and is balanced by strong 
capillary tensions associated with the curvature of the fluid-fluid interface. 
Consequently, expression (23) introduces a new term associated with the 
contact line Ft- 

Creeping flow of a Newtonian fluid is unrealistic at the corner of a moving 
contact line. Nevertheless, the system that is based on Eq. (4), boundary 
conditions (5), (6), adherence assumption on solid surfaces together with the 
dynamic Young-Dupre relation (23) for the apparent contact angle which 
accommodates the behaviour of the non-Newtonian fluid domain in the im- 
mediate vicinity of the contact line poses a problem of slow fluid motion. 
This result, in agreement with experiments and molecular investigations, is 
brought out by means of the analytic representation of the streamlines on 
Newtonian behaviour; however, other non-Newtonian fluid behaviour is used 
to arrive at bounded dissipative functions near the contact line. 

For slow movements we are able to give a model which provides answers 
to the previous questions: 

- The contact line is a non material line and acts in a similar way as a 
shock line. 

- The velocity fields are multivalued on the line. 

- The paradox of infinite viscous dissipation is removed. 

- Adherence and boundary conditions on surfaces and interfaces are pre- 
served, but a dynamic Young-Dupre relation derived by the virtual work 
principle yields the apparent dynamic contact angle as an implicit function 
of the contact line celerity. The apparent dynamic contact angle is the only 
pertinent Young angle from a continuum mechanics point of view. 

- For partial wetting, and for a sufficiently small capillary number, the 
concept of line friction is associated with the apparent contact angle. 
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The contact angle hysteresis phenomenon and the modelUng of exper- 
imentahy well-known results that express the dependence of the dynamic 
contact angle on the celerity of the line are important phenomena. In part 
2, Eq. (23) allows us to obtain an explanation of the contact-angle hysteresis 
in the advance and retreat of slowly moving fluids on a solid surface. 
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